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Abstract 

For a disordered system near the Anderson transition we show that the 
nearest-level-spacing distribution has the asymptotics P(s) oc exp(— As 2 " 1 ) 
for s ^> {s) = 1 which is universal and intermediate between the Gaussian 
asymptotics in a metal and the Poisson in an insulator. (Here the critical 
exponent < 7 < 1 and the numerical coefficient A depend only on the 
dimensionality d > 2) . It is obtained by mapping the energy level distribution 
to the Gibbs distribution for a classical one-dimensional gas with a pairwise 
interaction. The interaction, consistent with the universal asymptotics of the 
two-level correlation function found previously, is proved to be the power-law 
repulsion with the exponent —7. 
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It is well known [jTJ-Q that a statistical description of energy levels of quantum disordered 
systems in the metallic phase is provided by the random matrix theory (RMT) Its most 
important characteristic is the repulsion between the energy levels at any scale. In the 
Anderson insulator phase, the energy levels are independent and described by the Poisson 
statistics, provided that the appropriate states are separated by the length exceeding the 
localization length. 

It has been conjectured [|,|| that a universal statistical description is possible also in 
the critical region in the vicinity of the Anderson transition between the two phases. The 
dimensional scaling estimation made in Ref. [fj for the variance (N 2 ) — (N) 2 of the number 
of energy levels in a given energy interval, has suggested that it is proportional to (N) , thus 
being different from the Poisson statistics only by a certain number. A different statistical 
characteristic, the nearest-level-spacing distribution, has been conjectured in Ref. || on the 
basis of numerical simulations to be some universal 'hybrid' of the Poisson distribution for 
large level spacings and the Wigner surmise (see below) for small spacing. 

However, it has recently been analytically proved that the universal statistics, exactly 
applicable near the Anderson transition point (mobility edge), is entirely new and drastically 
different from both the RMT and the Poisson limit. The variance of the number of levels in 
the energy interval E S> A (centered at the Fermi energy) has been found as 

(N 2 )-(N) 2 = J^y = J(N)\ 0< 7 <1. (1) 

Here A is the mean level spacing, (...) denotes the ensemble averaging, the coefficient 
and the critical exponent 7 depend only on the dimensionality d > 2, and /3 is determined 
by the class of symmetry (/3 = 1, 2, or 4 for unitary, orthogonal, and symplectic ensembles, 
respectively Q). Eq. (1) is exact at the mobility edge in the limit 

L -> 00, (N) = const > 1. (2) 

In the same limit, the metallic phase is exactly described by the RMT and the insulator 
phase by the Poisson statistics. 
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Thus, the spectral rigidity does not disappear at the mobility edge (in contrast to the 
insulating phase where the levels are independent and the variance equals (N)), but it is 
considerably weaker than in the metallic phase (where the fluctuations are suppressed and 
the variance is proportional to ln(iV)). Then it is naturally to expect that the nearest-level- 
spacing distribution near the Anderson transition is also different from those both in the 
metal and insulating phases. 

Indeed, we will show in this Letter that the asymptotics of this distribution at the 
mobility edge is given by 

P(s) oc exp(-A d /?s 2 - 7 ), s = uj/A > 1, (3) 

where u is a distance between adjacent levels, and is some numerical factor depending 
only on the dimensionality d. This is drastically different from both the Poisson distribution, 
P(s) = exp(— s), and the exact Gaussian asymptotics in the metallic phase |4]] 

P(s) ~ exp {-^ 2 (3s 2 ) (4) 

Note that the famous Wigner surmise P(s) = {us/ 2) exp(— its 2 / 4) (for (3 = 1) describes this 
asymptotics only approximately [Q. 

Both the universal variance (1), and the asymptotics of the distribution (3) result from 
the exact asymptotics of the spectral density correlation function at the mobility edge 

R(uj) = ^lu(E)u(£'))-l = -c d l3- 1 \x-x'\- 2+1 , x = e/A, \x - x'\ > 1 (5) 

where v(s) is the exact density of states at the energy e, uq is the mean density of states, Cd 
is a positive number depending only on the dimensionality d > 2. The asymptotics @ has 
been obtained in Ref. |7j by calculating all the diagrams (with accuracy up to a numerical 
coefficient) which turned out to be possible after taking into account the analytical properties 
of the diffusion propagator and certain scaling relations at the mobility edge. 

To derive the announced result (|3|), we will use the effective "plasma model" as suggested 
by Dyson || . In such a model, the level distribution is mapped to the Gibbs distribution of 
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a classical one- dimensional gas of fictitious "particles" with a repulsive pairwise interaction 
f(\si — Ej\) in the presence of a confining potential V(e) 

V({e n }) = Z- 1 exp [-(3H({e n })} , (6) 

w(K» = E/(N-^l) + E^) (7) 

i<j i 

Here Z is the partition function and (3 plays a role of the inverse temperature. For 
f(\Ei — Ej\) = In |£j — , Eqs. (0) reproduce exactly the level distribution in the 
RMT, with /3 depending on the symmetry class as described after Eq. (1). The choice 
V(ei) = e\ for the confinement potential leads to the Gaussian ensembles but other choices 
are also possible |||| . In the metallic phase such a description is exact for the energy sepa- 
ration \si — Ej\ < E c = %/td where td = L 2 /D is the 'ergodic' time necessary for an electron 
to diffuse across the system. For t < td, i.e. \ei — £j\ > E c , the level statistics is completely 
different from that of the RMT. However, it is also described with the Gibbs distribution 
®j ©5 albeit with the pairwise interaction / having the power-law asymptotics [10]. 



At the mobility edge E c ~ A (g = E c /A is a dimensionless conductance). Therefore, 
the energy separation of a few A is already outside of the RMT region of validity. We will 
show that the asymptotics of the correlation function (|5|) at the mobility edge is described 
correctly by the Gibbs model with the power-law interaction 

f(\x-x'\) = ^^cot(7r 7 /2) |x-x'r 7 , 0<7<1, x = e/A. (8) 
2nc d 

(Naturally, this interaction is different from that in Ref. [Hj, where 7 > 1, which describes 



the nonuniversal level statistics in the metallic phase at the scale u 3> E c .) Before proving 
this, we will show how the form of the pairwise interaction governs the asymptotics of P(s). 

The distribution P(s) describes the probability to find the nearest adjacent level at the 
distance s = u/A from a given one. It is equivalent to the probability to find a "gap" of 
the width s (i.e. region that contains no "particles") in the Gibbs model. This probability 
is obtained from Eq. (|) as 

P(s) = exp {-P(F S - F )\ (9) 



where F s is the free energy of the one-dimensional gas (0) distributed along the straight 
line with the gap s around its center. For s>l, one introduces a continuous density p s (x) 
to describe such a distribution. Then, in the mean-field approximation (MFA) F s may be 
expressed as the functional 

F s = \l dx [ dx' p s (x)p s (x')f(\x-x'\) + ( dxp s (x)V(x), (10) 

l J\x\>^ J\x'\>% J\x\>§ 

where p s (x) obeys the mean- field (MF) equation 

[ dx' p s (x')f(\x-x'\) = -V(x) - p., \x\ > s/2. (11) 

Here p s arises from the "particles" number conservation (corresponding to the level number 
conservation in the original quantum disordered system), 



x\>i 



roc 

dx p s {x)dx = \ po(x)=Af. (12) 



Taking s = in Eqs. ([[(]) and (|TTD, one finds the density po(x) and free energy F for a 
homogeneous distribution. 

Equations similar to (|K]) and ( |TTD have been derived for a circular ensemble (a classical 
gas with the log \x — x'l" 1 interaction confined to a circular wire) by Dyson ]S| (see also Ref. 
Q) who has also found the corrections to the MF solution (allowing for the entropy term 



added to the functional (|T0l) and for the discreteness of the original distribution) lead to 
a linear in s contribution to the difference F s — Fq. For s ^> 1 this contribution is small 
compared to the leading quadratic term, Eq. (f|). We will consider only the terms leading 
in the s>l limit which are described by the MFA. 

For the circular ensemble there was no need in the confining potential V(x). We use 
the linear ensemble that is more convenient for relating the interaction /(|£j — £j\) to the 
correlation function (|^). Equation ( |TT| ) with the weakly singular kernel (||) can be solved 



for any V(x). This exact solution shows a strong dependence of both p s and po (and thus 
F s and F ) on V(x) (see Eq. ( |TTD for p ). However, it is easy to show that the difference 
F s — F , and thus the distribution (||), does not depend on V for s ^> 1 in the limit Af — > oo. 



Furthermore, the asymptotics of this universal distribution may be found, with accuracy up 
to a numerical coefficient, without knowing the explicit solution to Eq. (|TT|). 

The explicit dependence of F s — F on V(x), Eq. fllCf), is excluded straightforwardly with 
the help of Eqs. (|TTD, (|T2j). Then, after some transformations using the fact that the change 
in the "chemical potential" due to the gap formation fi s — /i ~ s/Af < 1, one finds with 
the accuracy up to s/M that 

F s - F Q = 

Iff 1 rs/2 rs/2 

dx / dx' Sp x Sp(x')f(\x — x'\) + - / dx I dx' po(x)po(x')f(\x—x'\) (13) 



2 J\x\>% J\x'\>2 2 J-s/2 J-s/2 

where Sp(x) = p s (x) — po(x) decreases rapidly for x ^> s. The function 5p(x) obeys the MF 
equation: 

rs/2 

dx' Sp(x') /(|x— x'\) = / dx'po(x') f(\x— x\), \x\ > s/2, (14) 

x'|>f J-s/2 



which follows from Eq. (11), if one neglects the small term p s — f^o ~ s/M . The homogeneous 
level density po(x) in Eqs. (|l"3|), (|T^ ) still depends on ^(x), Eq. (|Tl]). However, for \x\ < 
s/2 ^ M, this dependence is negligible, and in the limit M — ► oo one finds po = 1 ( m units of 
1/A). Now it is clearly seen from Eqs. flT3] ) and (|TJ) that the quantity F s —F = — In P(s) 
is determined by the interaction f(\x — x'\) only. 



Equations ([13]) and (|14]) are valid for an arbitrary long-ranged interaction x'\). For 
the case of the power-law interaction (H) (with < 7 < 1), one may rescale x — > sx and 
x' — > sx' to find that the solution of Eq. (|i~4T) (with p =const) has the form 5p(x) = <p(x/s), 
where (f(x) is a universal function independent of s. Substituting this solution to Eq. (p~3|) 
we arrive after the same rescaling at: 

F s -F = -/?- 1 lnP(s) = A 7 s 2 - 7 (15) 

where the universal constant depends only on the power of the interaction. Calculating 
this constant for the limiting case (7 = 0) of the logarithmic interaction in Eqs. (0), ( p^|) 



one reproduces the asymptotics (Q) known from the RMT || which leads to the announced 
result Eq. (|). 
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To prove that the interaction (H) reproduces the correlation function fl5|), we use the 
relationship ||: 

*>**~r l W)- (16) 

For the interaction / = I 7 and an arbitrary v(x) = V(x) + /io, the solution to Eq. 

(p|) with s = is found, using the methods described in Ref. |TlJ, as follows: 



7-1 



^ , cos 2 (7r 7 /2)(x + J D)^ , 
Po (x) = - — ^ 5 7, ^ x 



7T (Xry 



x T-l 1° dt{t + Df-\t-x)^ L f dTir + Dy-^it-r) 2 ^ v(r)\ , (17) 
ax [ at J-d J 

Here B is the Euler function, D is the band edge that may be found from Eq. fll2|) and 

tends to infinity when A/" — > oo. Taking the variational derivative ([16]), i.e. substituting 

— (3~ 1 5(t — x 1 ) for v(t) in Eq. (|T7D, one finds in the limit D — > oo: 

R(x, x') = -Z?" 1 ^^ cot \x - x'y~\ (18) 



27ra 7 V 2 / 

So the Gibbs model with the power-law interaction results in the asymptotics (|5|) of the 
correlation function. Comparing Eqs. (|5]) and (0), we obtain Eq. (]|). 

Note that for all the three universal statistics, in the metal and insulating phases and 
at the mobility edge, a simple relation holds between the variance of the level number 
fluctuations in the limit (|2]), and the asymptotics of the nearest-level-spacing distribution. 
Namely, if the variance proportional to (iV) 7 , then — lnP(s) oc s 2-7 . The linear in (N) 
variance is forbidden at the mobility edge [0 by the exact sum rule that is due to the 
conservation of the total number of states N ' . Therefore, the Poisson (i.e. linear in s) 
asymptotics of P(s) is equally forbidden at the mobility edge. Finally, following Ref. ||, we 
note that for s C 1 the distribution P(s) shows at the mobility edge the same behaviour as 
in the metallic phase, P(s) ~ s 13 , which follows from the general symmetry theorem proved 
by Dyson 0. Then, the whole distribution could be described by the following surmise: 

P{s) = Bs p exp {-Adps 2 -" 1 ) (19) 

where B is found from the normalization conditions. 
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